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Quantum spin tunneling and Kondo effect are two very different quantum phenomena that pro¬ 
duce the same effect on quantized spins, namely, the quenching of their magnetization. However, 
the nature of this quenching is very different so that quantum spin tunneling and Kondo effect 
compete with each other. Importantly, both quantum spin tunneling and Kondo effect produce 
very characteristic features in the spectral function that can be measured by means of single spin 
scanning tunneling spectroscopy and allows to probe the crossover from one regime to the other. We 
model this crossover, and the resulting changes in transport, using a non-perturbative treatment of 
a generalized Anderson model including magnetic anisotropy that leads to quantum spin tunneling. 
We predict that, at zero magnetic field, integer spins can feature a split-Kondo peak driven by 
quantum spin tunneling. 

PACS numbers: 


I. INTRODUCTION 

Quantum spin tunneling (QST) and Kondo effect are 
two ubiquitous and widely studied^. 3 phenomena in the 
broad field of nanoscale magnetism. They both turn a 
spin system with a doubly degenerate ground state into a 
system with a unique ground state with null magnetiza¬ 
tion. QST affects quantized integer spins with magnetic 
anisotropy, such as single molecule magnets^, magnetic 
impurities in insulators^ and magnetic adatoms^ and 
molecules^ on surfaces. The Kondo effect is most of¬ 
ten associated with half-integer spins, but it has been 
observed in a variety of integer spin systems, such as 
quantum dots with an even number of electrons^, var¬ 
ious integer spin magnetic molecules^— and molecular 
oxygen (spin S = 1) adsorbed on goldii. 

The Kondo effect arises when a local spin is exchange 
coupled to itinerant electrons that respond dynamically 
to screen the magnetic moment of the impurity^ This dy¬ 
namical response leads to a resonance in the local den¬ 
sity of states at zero energy, that emerges as a zero bias 
Fano feature in the transport spectroscopy curves, G(V), 
where G = The observation of the Kondo reso¬ 

nance in individual magnetic atoms^ and molecules^ by 
means of Scanning Tunneling Microscope (STM) spec¬ 
troscopy has been reported numerous times in the last 
two decades. The natural energy scale that character¬ 
izes the Kondo effect is roughly given by the width of 
this resonance which depends on the tunneling rate T for 
electrons between the localized atomic orbitals and the 
extended states of the surface, and the charging energy 
of the atom, Un = E(N + 1) — E(N). 

QST can occur for integer spins with negative domi¬ 
nant uniaxial anisotropy D < 0: In this case the ground 
state of a spin S is doubly degenerate and consists of 
the two states with opposite and maximal spin projec¬ 
tion m z = ±S, which are separated by an energy barrier 
~ DS 2 . A finite in-plane magnetic anisotropy E then al¬ 


lows quantum tunneling between the two spin states, lift¬ 
ing the degenerac y 13 ’ 14 of the ground state by the tunnel 
splitting Ao- Thus the ground state is a linear combina¬ 
tion of two spin states with opposite spin projection and 
null magnetizatio n 15 ’ 16 . The dynamical spin response 
function acquires a pole at hui = A 0 so that QST can be 
probed spectroscopically. Whereas in large spin molecu¬ 
lar magnets Ao is so small that it can only be inferred 
indirectly^, for systems with S = 1,2, such as mag¬ 
netic adatoms and small magnetic molecules deposited on 
conducting substrates, the spin excitations, and thereby 
the QST splittings, have been resolved directly by means 
of inelastic electron tunneling spectroscopy (IETS) us¬ 
ing STM^. The corresponding G{V) spectra show step 
features and, in general, no Kondo peak. It is cus¬ 
tomarily assumed that these excitation energies are a 
property of the atomic/molecular quantized spin, weakly 
dressed^i£r-i£ by its Kondo exchange coupling. 

Whether a given magnetic atom or molecule will show 
a stepwise G(V ) spectrum^ - — or a zero bias Kondo reso¬ 
nance depends on the strength of the Kondo exchange 
J which is controlled by the ratio T/U. Thus, FePc 
molecules, that in vacuum have S = 1, display a Kondo 
feature when deposited on Au(lll)^ and inelastic steps 
when deposited on oxidized Cu(110)£. Moreover, the 
joint observation of a zero bias Kondo resonance together 
with stepwise inelastic spin excitations has been reported 
for individual magnetic atom o 17 ' 19 . Importantly, it is 
possible to devise experiments^ ’ 17 ’ 20 in which the Kondo 
interaction could be tuned, making it relevant to address 
the question of how the G(V) spectra evolve from the 
weak to the strong coupling regime. The magnitude of 
the quantum spin tunneling splitting can also be modu¬ 
lated by application of a magnetic field along the hard 
axis directio n 13 ’ 14 . In this work we address how the com¬ 
petition between Kondo screening and QST affect the 
STM inelastic conductance and we predict a new physi¬ 
cal phenomenon, the splitting of the Kondo peak at zero 
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FIG. 1: (a) Schematic model of experimental setup for mea¬ 

suring excitation spectra of a magnetic adatom on a surface 
with STM tip. (b) Schematic energy spectrum of adatom GS 
multiplet for S = 1 for negative uniaxial anisotropy D < 0 and 
finite in-plane anisotropy E. (c) Same as (b) but for 5 = 2. 
(d) Comparison of PP spectrum A m (ca) = A m {oj)/ uj) for 
the four lowest excitations (bottom) and the real electron 
spectrum (top) for 5 = 2, D = — 1.55meF and E = 0.35meV. 
The vertical grey lines indicate the positions of the PP ener¬ 
gies E^ w.r.t. the GS PP energy Eg. 

magnetic field due to quantum spin tunneling. 

Both the inelastic steps^i and the Kondo features^ - — 
can be described using a Kondo Hamiltonian where the 
atomic spin is described with a single-ion quantized 
spin interacting, via exchange, with the conduction elec¬ 
trons of the surface. In particular, the interplay be¬ 
tween magnetic anisotropy and Kondo screening has been 
thoroughly studied using the Kondo mode l 17 ’ 18 i 22 v 21. 
Here we present a more general approach, based on a 
non-perturbative treatment of a multi-orbital Anderson 
model for the adatom coupled to the substrate, that per¬ 
mits to include atomic charge fluctuations that are effec¬ 
tively frozen in the Kondo model. As we show below, 
valence fluctuations have a similar effect on the spec¬ 
tra as Kondo exchange. In addition, density functional 
theory (DFT) calculations show that often charge is not 
quantized in magnetic adatom system s 25 ’ 26 . 

II. MODEL AND METHOD 

We consider tunneling between an STM tip and a mag¬ 
netic adatom coupled to a surface as shown in Fig. [Qa). 
Assuming weak coupling to the STM tip (tunneling 


regime) 27 the low-bias conductance can be directly re¬ 
lated to the adatom density of states (technically, the 
many-body spectral function) p a (u>), 

S(V) = ^-J2 T a P P a (eV) (1) 

a 

where r^ lp = is the (weak) tunneling rate of 

electrons between the adatom orbitals a and the STM 
tip. We have assumed the DOS of the STM tip p t i p is 
energy independent around the Fermi level. Note how¬ 
ever, that in general different orbitals couple differently 
to the STM tip so that the contribution of the individual 
channels to the total conductance may differ. We neglect 
direct tunneling into surface states in eq. dTJ. This is a 
good approximation when the magnetic atoms are sep¬ 
arated from the metallic surface by a decoupling insu¬ 
lating layer, such as Cu 2 N/Cu(100)liiid^, CuO/Cr£ and 
h-BN/Rh(lll)^. This approximation does not capture 
the Fano interference effect relevant^ when the tip-atom 
channel interferes with the direct tip-surface tunneling 
path^S. 

We describe the magnetic atom on the surface by a 
multi-orbital Anderson model, 

E = 'Himp + Hbath + Vhyb (2) 

where the Hamiltonian of the Anderson impurity site 
T~iim P describes the strongly interacting 3d-levels that 
yield the spin of the magnetic atom, and includes a term 
that accounts for magnetic anisotropy: 

%imp = -f- ^ ^ U TIolo Tla.' o' ^ S a * S a f 

olo^ol' a' ol^ol' 

+ DSl+E{S 2 x -S 2 v ) (3) 

ed are the single-particle energies of the d-levels, Nd = 

Sa a 1S the number operator for all d-levels a = 

1, h at7 = <F a<r d a(r is the number operator of 

an individual d-level a with spin a , U is the effective 
Coulomb repulsion, Jh the Hund’s coupling, and S a 
measures the total spin of an individual d-level a, i.e. 
S a = Ylaa' datjTaa'daa'■ The crystal field splitting of 
the d-levels together with the spin-orbit coupling (SOC) 
gives rise^ to magnetic anisotropy (MA) which in our sim¬ 
plified model is taken into account by the effective spin 
Hamiltonian given by the last term of ([3]) where D is the 
uniaxial anisotropy and E the in-plane anisotropy^. 

The second term in m describes the conduction elec¬ 
tron bath in the surface: 

*^bath = E £ kac\ aa Ckaa (4) 

fe,a,cr 

The third term in ([2]) is the so-called hybridization term 
which describes the coupling between the impurity and 
the conduction electron bath: 

Vhyb = 'y ' Vka(c\. aa d aa + d^Ckao) (5) 

fc,a,cr 
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Integrating out the bath degrees of freedom one obtains 
the so-called hybridization function: 


=e 

k 


\v ka \ 2 

U> + H - £ ka + IT] 


( 6 ) 


Its (negative) imaginary part r a (w) = —ImA byb (w) de¬ 
scribes the single-particle broadening of individual impu¬ 
rity levels a due to the coupling to the conduction elec¬ 
trons. Note that we have assumed here that each impu¬ 
rity level a couples to different conduction electron states 
(labeled by a) in the substrate so that each impurity level 
has its own bath. This assumption is justified because of 
the different symmetries of the d-orbitals. In the case of 
coupling of two orbitals to the same conduction electron 
states, off-diagonal elements in the hybridization func¬ 
tion Aj a , (ui) would occur, describing substrate mediated 
hopping between impurity levels. In adatom-substrate 
systems these off-diagonal elements are often either zero 
or very small^A 

We now solve the Anderson model (2) within the so- 
called One-Crossing Approximation (OCA)22r—. The 
first step is an exact diagonalization of the isolated im¬ 
purity Hamiltonian ©: 

'Himp = E E ™ l m ) ( m l ( 7 ) 


The many-body eigenstates \m) are eigenstates of the to¬ 
tal number of electrons of the impurity, i.e. N& |to) = 
N m |to), and the total spin 5 2 of the impurity, i.e. 
S 2 \m) = S m (S m + 1) |to). For a ground state (GS) with 
integer spin 5, and magnetic anisotropy with negative 
uniaxial anisotropy D ^ 0 and finite in-plane anisotropy 
E ^ 0, the (25 + 1) degeneracy of the GS multiplet is 
completely lifted. This is schematically shown on the 
right hand side of Fig. [T| for 5 = 1 and 5 = 2 and 
D < For D < 0 and E = 0 the GS is doubly de¬ 
generate with the GS doublet having the maximal spin 
projection m z = ±5. A finite in-plane anisotropy E al¬ 
lows for quantum tunneling between both spin directions 
thus lifting the degeneracy of the GS doublet, which now 
becomes split by Ao, the bare quantum spin tunneling. 
The quantum states of the split doublet are thus linear 
combinations |±) ~ \m z = +5) ± | m z = —5). 

In the next step a diagrammatic expansion of the 
many-body eigenstates | to) of the (isolated) impurity 
Himp in terms of the hybridization © is developed. Vh y b 
connects eigenstates | to) and | n) of Hi mp with occupa¬ 
tion numbers differing by one (N m = N n ± 1, see (IA1I) 
in App. It is these fluctuations between the impurity 
GS and excited states with one more or one less electron 
that give rise to the Kondo effect. To this end one intro¬ 
duces so-called pseudo-particles (PPs) m corresponding 
to the many-body eigenstates |to). The full propagator 
of such a PP to can be written as 


G m (w) 


1 

uj A E m S m (cu) 


( 8 ) 


where £ m (w) is the PP self-energy which describes the 
renormalization (real part) and broadening (imaginary 
part) of the PP m due to the interaction with other PPs 
to' mediated by the conduction electron bath [see eq. 
mi in A PP- S- —A is the chemical potential for the 
PPs which has to be adjusted such that the total PP 
charge is constrained to one [see eqs. (1A4IA6I) in Apd.IA], 

OCA consists in a diagrammatic expansion of the PP 
self-energies £ m in terms of the hybridization function 
A byb to infinite order but summing only a subset of di¬ 
agrams (only those involving conduction electron lines 
crossing at most once). This leads to a set of cou¬ 
pled integral equations for the PP propagators and self¬ 
energies that have to be solved self-consistently. Once 
the OCA equations are solved the real electron spectral 
function p a {oj) for the impurity levels entering equation 
© for calculating the conductance spectrum is obtained 
from convolutions of the PP spectral functions A m (uj ) = 
—ImG m (w)/7r [see eqs. (lAlfillAlTl) in App. [A] which fea¬ 
ture sharp resonances at the renormalized many-body 
energies = E m + Re ’S rn (E* n ). The differences be¬ 
tween the renormalized energies of the excited states E* n 
and the GS Eq yield the real electronic excitations as 
shown Fig. 1(d). 

OCA captures both the weak and strong coupling 
Kondo regimes, and has shown to produce reliable spec¬ 
tra for the single-orbital Anderson model, as long as the 
temperatures are not too low compared to the Kondo 
temperature^. For the general multi-orbital situation 
considered here, benchmarking of spectra is difficult since 
Numerical Renormalization Group^l is computationally 
too demanding to be applied. Using Continuous-Time 
Quantum Monte-Carlo, it has been showed that OCA 
is markedly superior to the simpler Non-Crossing Ap¬ 
proximation (NCA) for dealing with multi-orbital An¬ 
derson models, although certain sum rule violations 
are found. Very importantly, OCA has shown excel¬ 
lent agreement with experiments in very complex rnulti- 
orbital system o 17 i 39 i 40 . More details about the OCA 
method are given in App. [© 


III. RESULTS 

We assume that the spin 5 of the impurity is carried 
by 25 active impurity levels close to half-filling. For all 
impurity levels we assume the same energy-independent 
broadening T. The important energy scale is the Kondo 
exchange coupling J ~ T/U. which can be varied by either 
changing T or the charging energy IA which depends on 
the Coulomb repulsion U. Here we choose to vary T 
and fix the Coulomb repulsion to U = 4eV. The Hund’s 
coupling whose main effect is to favor the formation of a 
high-spin ground state is fixed to Jr = leV. 

We first study the competition between Kondo quench¬ 
ing and QST for a magnetic atom with spin 5 = 1, at 
the point of electron-hole (e-h) symmetry. In this case 
the competition is controlled by two energy scales, the 
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FIG. 2: Results for spin S = 1 and for D < 0, and tempera¬ 
ture kT = O.lmeV. (a-d) Spectra pd{oj) in dependence of dif¬ 
ferent system parameters: (a) Changing E while keeping fixed 
Ai = \D\ + E =10meV and F/n = lOOmeV. (b) Changing F 
for D = —9.9meV and E = O.lmeV. (c) Temperature depen¬ 
dence for r/7T = lOOmeV, D = — 9.75meV and E = 0.25meV. 

(d) Same as (c) but for D = — 9.9meV and E = O.lmeV. 

(e) Effective QST Ao as a function of T. (f) Spectra for 
D = — 9.75meV, E = 0.25meV and r/7r = lOOmeV for differ¬ 
ent energy level shifts Std- In parenthesis the total occupancy 
Nd of the impurity is given. 


in-plane anisotropy E, that drives the QST, and T, that 
favors Kondo coupling. We assume D < 0 so that, for 
E = 0, the ground state is the doublet of eigenstates of 
S z with m z = ±1 and the state with m z = 0 is the first 
excited state [see Fig. QJb)]. Because of the two-fold de¬ 
generacy of the GS the Kondo effect can take place [top 
curve in Fig. [21(a)]. Spin-flip events for the GS doublet 
to = ±1 occur via the excited state m — 0 and hence are 
reduced by a factor of J/D compared to the Kondo ex¬ 
change J of the corresponding spin-1 Kondo model with¬ 
out anisotropy. Note that such a spin-flip process involves 
the simultaneous scattering of two conduction electrons. 
This is possible since we are considering a multichannel 
situation as each impurity orbital is connected to its own 
bath. In the situation of just a single screening channel, 
the spin-flip between the m = ±1 GS doublet would be 
inhibited when D becomes bigger than the Kondo tem¬ 
perature, leading to a split-Kondo featur o 41 ’ 42 . 

The in-plane anisotropy term E > 0 produces the 
QST that leads to a splitting of the to = ±1 doublet, 
so that the bare excitation energies are Ao = 2 E and 
Ai = \D\ + E. At weak coupling (small T), the spec¬ 


tra show two steps corresponding to inelastic spin transi¬ 
tions between the renormalized spin levels [bottom curve 
in Fig. Ufa)]. Our calculations show that the effect of 
decreasing E [Fig. [21(a)] and increasing T [Fig. [21(b)] is 
similar. In both instances the renormalized QST split¬ 
ting Ao decreases, and the line shapes evolve from square 
steps at small T or large E to a characteristic triangular 
shape, very often seen in experiments!!!!, that can only 
be captured in part when going beyond second order per¬ 
turbation theory in the Kondo exchang o 23 ’ 24 . Whereas 
the reduction of the renormalized QST splitting Ao as 
E decreases is trivially accounted for by the fact that 
Ao = 2 E, the red-shift renormalization of the spin ex¬ 
citation energies due to Kondo coupling [Fig. [21(e)] is a 
many-body effect, in line with previous results!! -18 ’ 29 ’ 30 . 
In our approach the renormalization of the spin excita¬ 
tions ultimately originates in the renormalization of the 
many-body energies E^ by the real part of the PP self¬ 
energy S m (w), see eq. ©. We would like to stress at 
this point that the splitting of the Kondo peak is exactly 
given by the effective QST A 0 . Hence Fig. [21 e) really is 
a prediction for the splitting of the Kondo peak by QST 
which can be measured experimentally. 

At intermediate coupling, the renormalization of the 
line shape and energy of the lowest energy spin excitation 
turns them into a split Kondo peak. This is one of the 
important results of this work: In the absence of a mag¬ 
netic field the competition between QST and Kondo ef¬ 
fect yields a split-Kondo state whose characteristic signa¬ 
ture can be measured by STM spectroscopy. As shown in 
Figs. [2Kc,d), the temperature dependence of the spectra 
in the case of intermediate coupling (split-Kondo) is quite 
different from that of weak coupling (step-like). For weak 
coupling [Fig. [21(c)], at low temperatures two clear steps 
are obtained, with slight triangular departures from the 
step-wise behavior, due to Kondo interactions, that are 
captured as well by perturbative calculations!!. As the 
temperature is increased, the low energy step is smeared 
out and, eventually, is no longer resolved. On the other 
hand, at intermediate coupling [Fig. [21(d)], a zero-field 
split Kondo peak is obtained at low temperatures. As 
we increase fcgT the QST splitting of the Kondo peak 
disappears due to thermal smearing, resulting in a single 
peak. At strong coupling, the effective QST Ao van¬ 
ishes [Fig. [21(e)], so that the ground state and the first 
excited state become effectively degenerate: Kondo cou¬ 
pling quenches QST splitting!!. At this point, the spec¬ 
tra show a Kondo peak [top curves in Figs. [2)(a,b)]. 

We now consider the effect of valence fluctuations. 
This effect is particularly important since DFT calcu¬ 
lations show that often the average occupation of the 
d shell is not quantize d 25 ’ 26 . In our model, we control 
the valence mixing by shifting the impurity levels by an 
amount Sej, taking the system out of the e-h symme¬ 
try point, leading to deviations from integer occupation 
number of the impurity shell. Valence fluctuations have 
a similar effect on the spectra as increasing the Kondo 
screening, leading e.g. to a considerable enhancement of 
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FIG. 3: Results for spin S = 2 with D = — 1.55meV and 
E = 0.35meV, and temperature kT = 0.2meV. (a) Spectral 
function Pd(u) for different values of single-particle broaden¬ 
ing T. (b) Temperature dependence of spectral function for 
r/7r = lOOmeV. (c) QST Ao as a function of single-particle 
broadening T. (d) Effect of charge fluctuations: spectra for 
r/7T = lOOmeV for different energy level shifts 5td- In paren¬ 
thesis the total occupancy Nd of the impurity levels is given. 


the Kondo peah^. Similarly, also the renormalization of 
the excitation energies and especially of the QST by the 
Kondo is enhanced by the introduction of charge fluctua¬ 
tions, as can be seen in Fig. [3]T). Increasing 5ed leads to 
the step features associated with spin excitations mov¬ 
ing to lower energies. Also the up-bending of the low¬ 
est energy steps associated with QST and the conversion 
into a split-Kondo peak is induced by increasing valence 
fluctuations. Different from the case when T increases, 
the spectra also become somewhat asymmetric as the e-h 
symmetry is broken by the valence fluctuations. 

The competition between Kondo effect and QST is also 
present for higher integer spin systems. In Fig. [3] we 
show results for spin 5 = 2, and magnetic anisotropy 
parameters D = — 1.55meV and E = 0.35meV chosen 
to reproduce spectra for Fe on Cu 2 N^. The energy level 
diagram, and the spin composition of the corresponding 
states, are shown in Fig. [TJ Here the QST splits both the 
m = ±2 ground state doublet, by an amount Ao oc yy 
as well as the m = ±1 doublet by 2 E. Thus, the mag¬ 
netic anisotropy completely lifts the degeneracy of the GS 
quintuplet leading to step features in the calculated spec¬ 
tral functions [Fig. [3ja)] that, for small values of T, re¬ 
semble those measured for Fe on Cu 2 N^, and correspond 
to inelastic spin excitations from the ground state. As 
r is increased, the inelastic spin excitation step features 
both move to lower energies and become broader, indicat¬ 
ing both the renormalization of the associated excitation 
energies and the decrease of the lifetimes due to the ex¬ 
change coupling. As in the 5 = 1 case, the step features 
associated with QST develop a triangular shape with in¬ 
creasing r and finally turn into a split-Kondo peak. As 
r is increased, the effective QST Aq is renormalized but 


remains finite in the range of physically reasonable val¬ 
ues of T (up to leV) considered here [see Fig. [Htc)] . We 
note that generally the Kondo coupling is also weaker for 
5 = 2 than for 5 = 1 since higher order processes are nec¬ 
essary to screen the spin. For a somewhat smaller E (not 
shown) the splitting of the Kondo peak can also vanish 
for spin 5 = 2. Fig. [3[b) shows the temperature de¬ 
pendence of the spectrum for the case of relatively small 
coupling r. As the temperature is lowered, initially flat 
steps appear in the spectrum, which become increasingly 
triangular as the temperature is lowered further. As in 
the case of spin 5 = 1, increasing the valence fluctua¬ 
tions by detuning the system from ph symmetry has a 
similar effect as increasing T [see Fig. [31(d)] , leading to 
a red-shift of the spin excitation energies. The detuning 
from ph symmetry also leads to an up-bending of the ini¬ 
tially flat step features associated with QST and to the 
conversion into a split-Kondo peak. 

The picture that emerges from our calculations is the 
following. An atom with integer spin 5, described with 
Hamiltonian ([3]) , is a closed quantum system whose quan¬ 
tum ground state has a built-in coherence between the 
two classical ground states, with S z = ±5, associated 
to the QST splitting Ao, that would determine the fre¬ 
quency of the Rabi flops of the magnetization if the atom 
was initially prepared in an eigenstate of S z . As the cou¬ 
pling to the surface electrons is turned on, the atomic 
spin behaves like a quantum open system. This results 
both in the renormalization (reduction) of this Rabi fre¬ 
quency as well as spin relaxation, leading to a broaden¬ 
ing of the steps. For sufficiently large coupling the QST 
splitting can be completely quenched, in line with previ¬ 
ous results^, and quantum coherence between the states 
with opposite S z is lost. At that point Kondo screening 
takes over, and the zero bias Kondo feature appears. 

Interestingly, the evolution between the strong cou¬ 
pling Kondo regime and the weak coupling with step-wise 
excitations is continuous. At intermediate couplings the 
spectral functions show a split Kondo peak, that resem¬ 
bles the Zeeman split Kondo peak, but is driven by QST 
instead. The strength of the Kondo coupling is controlled 
both by r and by the departure from the e-h symmetry 
point, that we change by tuning Sed■ Importantly, both 
the weak coupling picture with step-wise excitations as¬ 
sociated with spin transitions, and the strong coupling 
Kondo regime, are also obtained in the case when the 
charge on the atom is not quantized, as suggested by 
DFT calculation s 25 : 26 . 


IV. CONCLUSIONS 

In summary, we have studied the competition between 
two important physical phenomena that affect integer 
spins, namely QST and Kondo screening. Our calcula¬ 
tions permit to trace the evolution from the weak Kondo 
coupling regime, where the stepwise dl/dV spectra are 
renormalized, resulting in the shift and the broaden- 
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ing of the spin excitation energies, to the strong cou¬ 
pling regime, where the zero bias Kondo peak appears. 
This accounts for several experimental observation ;) 17 i 20 . 
For strong Kondo coupling, QST can be completely 
quenched. Importantly, for the intermediate coupling 
regime we predict a new physical effect: for B = 0 an 
energy split Kondo peak can arise, because of quantum 
spin tunneling splitting, in analogy with the B ^ 0 Zee- 
man split Kondo peak, recently investigated more closely 
by spin-polarized STM.— 
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where |0) is the PP vacuum. The auxiliary fields a 
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Appendix A: One-Crossing Approximation 

The starting point are the eigenstates | to) and cor¬ 
responding eigenenergies E m of the isolated impurity 
Hamiltonian dTJ). In terms of the impurity eigenstates 
| m) we can rewrite the hybridization term Vhyb as 

Vhyb = EE Ffca (|m.) (m| d) a(7 | ti) (n| Ck a(T T h.c.) 

m,n k,a,a 

(Al) 

Vhyb connects eigenstates of Hi m p with different occupa¬ 
tion numbers, i.e. N m = N n ±1. It is the fluctuations 
between the GS and excited states of the atom give rise to 
the Kondo effect. Since the electron-electron interaction 
on the impurity (U and Jh ) is generally large compared 
to the hybridization (I4 a ), a perturbative treatment in 
terms of the latter is justified. 

We are now going to develop a perturbation theory for 
the many-body eigenstates | m) of 'Himp in terms of the 
hybridization term Vhyb- In order to proceed we asso¬ 
ciate so-called pseudo-particle (PP) field operators a^, 
a m with the many-body eigenstates | to) of Himp: 

|m) = ajn |0) and a m \m) = |0) (A2) 


Q — ^ \ — 1 (A4) 

m 

The PPs are related to the real electrons by: 

d a a = ^2 ( m\ d aa | n) (A5) 

m,n 


In the PP picture the Hamiltonian of the isolated atom 
becomes 



where the last term is a Lagrangian constraint impos¬ 
ing the afore mentioned conservation of the PP charge 
Q. The Lagrange multiplier A can be seen as a (nega¬ 
tive) chemical potential for the PPs. Imposing the con¬ 
straint via the Lagrange multiplier is done e.g. in the 
Slave-Boson Mean-Field Approximation of the Anderson 
model.— However, for the diagrammatic expansion it is 
more convenient to work in the grand-canonical ensem¬ 
ble with respect to the PP charge Q. The expectation 
values of a physical observable A can be calculated in the 
grand canonical PP ensemble and then projected to the 
physical subspace (Q = 1) via the Abrikosov tricks^ 


(A) 


lim 

A->oo (Q) x 


(A7) 


In the PP picture the hybridization term (IA1I) be¬ 
comes: 


Vhyb — y ] y ^ Vfc a {^D m n al^anCkaa + -Orrm C feacr a m a ii) 

m,n k,a,cr 

(A8) 

where we have introduced the matrix elements D^ n = 
(to| d atT \n) and = (m| d' arj |n). The hybridiza¬ 

tion term Vhyb is now treated as a perturbation to 
the Hamiltonian of the uncoupled impurity and bath 
Ho = Himp + Hbath ■ The introduction of PPs obeying 
(anti-) commutation relations allows us to make use of 
the machinery of quantum field theory to develop a dia¬ 
grammatic perturbation expansion since Wick’s theorem 
applies. 
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We now introduce the bare PP propagators G^(r) = 
— (T T a m (r)a] Tl ( 0 ))o which will be denoted by dashed 
lines: 




1 

iuj — A — E m 



(A9) 


The bath electron propagator is denoted by a full line: 

Qkaa (iw) — 7 = - > (AlO) 

e fea kaa 


Expanding the full PP propagator G m (r) = 
(T r a m (r)a] n (0)) in terms of the hybridization Vhyb, we 
can integrate out the bath electrons, and end up with an 
effective retarded two-particle interaction between PPs 
mediated by the bath electrons: 


Vhyb(r-r')= £ D^D%7 m ,x (All) 

a,a, 171 , 71 , 171 ' ,n' 


x a^(T)a„(r) 


^2 \Vka\ 2 gkaaij ~ t') 


L k 




The term in square brackets is called the hybridization 
function A a (r — r') which in the real frequency domain 
is given by eq. ©. 

The full PP propagators are denoted by double dashed 
lines. In terms of a PP self-energy E m (w) which captures 
the interaction with other PPs via Vhyb, the full propa¬ 
gator can be written as 


m iuj - A — E m - H m (iu) m 

(A12) 

In lowest order (2nd order in T4 a or 1st order A a ) the 
PP self-energy is given by the following diagrams: 


where r a (u;) = — ImA a (u;) is the single-particle broad¬ 
ening of the impurity levels a by the coupling to the 
conduction electron bath, and f{y) is the Fermi func¬ 
tion. Hence the NCA self-energy for a PP m is given by 
a convolution of the imaginary part of the hybridization 
function r a with the propagators of all other PPs m! that 
m is interacting with via the conduction electron bath. 

The OCA diagrams are second order in A a (u>) and 
involve crossing conduction electron lines: 



(A15) 


The algebraic expressions for the OCA self-energy are 
much more complicated than the NCA ones, and involve 
double convolutions of two hybridization functions with 
three PP propagators. The exact expressions can be 
found in the literature ! 34 ! 35 

The NCA/OCA equations are a set of coupled integral 
equations that have to be solved self-consistently since 
the self-energy of a PP m depends on the full propaga¬ 
tors of other PPs m!. Once the NCA/OCA equations 
are solved, the real impurity electron propagator can be 
determined by making use of the relation (IA5I) . Hence 
the real electron propagator can be calculated from a 
two-particle correlation function for PPs: 


ai7 aa 



m! m! 


(A13) 


The conduction electron line (full lines) aa correspond 
to the hybridization function A a (u;). 

The first diagram (“backward” diagram) corresponds 
to adding an electron to the impurity site, i.e. N m > = 
N m + 1 and the second diagram (“forward” diagram) to 
removing an electron, i.e. N m n = N m — 1 . The non¬ 
crossing approximation (NCA) consists in an infinite re¬ 
summation of these lowest order diagrams where conduc¬ 
tion electron lines do not cross (hence the name). Replac¬ 
ing the bare propagators by full propagators in the above 
self-energy diagrams, one obtains the NCA self-energy di¬ 
agrams: 


at7 aa 



vn! w! 


- E / ^r[\ D Zn '\ 2 /M r.M + 

m' aa 

+ \DZ’m\ 2 /(-*') r «M G m .{cc - v)\ (A14) 


G aff (r) = -(T T d a<J (T)dl a ( 0)) 

= - E 

171 , 71 , 171 ' ,n' 

(A16) 


In NCA the vertex correction A is neglected. It is then 
found that the real electron spectral function p ao .(w) = 
—TmG aa (oj)/'K can be calculated from a convolution of 
PP spectral functions A m (w) = -ImG m (u)/ir: 


p a M 


A r £/*.-*tL + .-*ix 

' ' mm' 

X l^ram'| 2 ^m( £ )^m'(w+£) (A17) 


where Q is the PP charge which can be calculated directly 
from the PP spectral functions: 

(Q) = f duje~ 0UJ ^2 A m (w) (A18) 

m 

Again the corresponding expression for calculating the 
real electron spectral function within OCA is much more 
complicated as it involves double convolutions of PP cor¬ 
relation functions and the hybridization function ! 34 ! 35 
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